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ABSTRACT 

In this paper a unified presentation is made on the results of var- 
ious investigators on the properties of random communication networks. 

These results are interpreted in such a way that the properties may be 
determined by using a digital computer with the application of the 
Monte Carlo method. The computer program is written and tested. Re- 
sults for some networks are compared with theoretical values. 
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I. RANDOM COMMUNICATION NETWORK 



In this paper the results of various investigators on the properties 
of random communication networks are reviewed. From these results a 
mathematical model is represented for a general random communication 
network . 

A. THE CONCEPT OF A RANDOM COMMUNICATION NETWORK 

A communication network is an aggregate of message centers that at- 
tempt to transfer information to one another over a wide variety of 
channels. The message centers could be fixed as the headquarters of a 
regional military command, mobile as ships, aircrafts, or satellites. 

The channels could be a radio link, a telephone line, or highway. In 
normal conditions, the topology of the network is known with certainty: 
the positions of ships at sea, the links between the headquarters with 
different units. Each center knows its relative position with respect 
to the other centers exactly. If a message is sent from a central 
command, the authorities know exactly how many and which units will re- 
ceive it. The network, in this case, is called a deterministic 
communication network. But during the hostilities, some units could be 
destroyed by the enemy without even having the chance to send the last 
word, some radio link could be jammed by the enemy, some telephone lines 
or some roads could be damaged without having any means to determine the 
extent of the destruction. Back at the headquarters, the authorities 
are faced with the uncertainty of the existence of the different units 
and channels, the effectiveness of an order sent to the battlefield is 
not as high as before. It is only probable that such unit gets the 
message. The communication network becomes a highly probabilistic, 
random netwrok in the event of a war. 
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One finds the same kind of randomness in many other branches of 
research; for example, the contagion of a disease among a population, 
an hereditary defect in a species, the spread of a stimulus in the 
nervous system, or the spread of information in society, to name a few. 
All these phenomena are very similar to the random communication net- 
work in that they have: 

a. an originator which originates a specific message 

b. several recipients which convey the received 'message' to the 
others via a medium (atmosphere, voice, physical contact, etc.) 

c. the uncertainty of contact between two individuals. 

All these phenomena could be then studied under the same topic of 
the random communication network. Researches in these phenomena give 
much insight in the problem of propagation message through a random 
network. 

As it is said, the communication between two message centers A and 
B could be assumed only with a certain probability. The latter measure 
is a function of many parameters which are for example, the time t 
when the message is sent, the time T since the message is received at 
station A, the distance d between stations, the electronic equipment of 
each station, to cite a few. Therefore one could write, in general: 

Prob (contact A to B) = p (A, B, t ,X. , d, . . .). 

Due to the complexity of the problem, without reasonable assumptions 
about the parameters to simplify the problem, one could not solve it. 
Following is the brief survey of the contributions of various workers. 
Each of them has solved partially the problem with different assumptions 
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about the parameters which determine the probability of contact between 
two stations, two individuals, whatever is the case. 

B. THEORETICAL DEVELOPMENT 

Bailey [1] studied the spread of an epidemic in an homogeneous 
population with size n where the probability of one new infection tak- 
ing place in the interval dt is given by 

Prob (one infection in dt) = y (n-y+l)dt, 
where y = the number of susceptibles . 

From the above equation, the mean m(t) of the number of infected people 
and the epidemic curve, (i.e., the rate of change with respect to time 

of the mean of the number of infected people) dm(t) are obtained to be 

dt 

m(t) = e- 10t (810t - 234) + e -18 *(4410 t - 902) + . . . 

z = dm(t)_ = e -l0tj- 8100 t _ 3156) + e - 18t (79,380t - 20650) + . . . 

dt 

Landau and Rapoport [2] assumed the probability of conta ct between 
any two individuals is the same for any pair and the probability of 
transmission of the disease depends both on the time t of the whole 
process and also onT, the time since the particular affected indivi- 
dual acquired the disease. Thus, 

Prob (one infection) = p(t,T) 

The differential equations giving the number of individuals who become 
affected and the rate are : 
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dx(t) = <*(N - x ) * fl p(t,t) 
dt L 




d z (t ) = _1 dx (t) = c< N ( 1 - z 
dt N d t 



)f *Q = p(t,t) ( t d_z p (t ,-t-A) 



N Jo d 



( 1 . 2 ) 



where N = number of individuals in the population 



x(t) = total of individuals who have become affected up to time t 
X 0 = number of individuals who become affected at t = o, the 
initial time of the process 

=frequency of contacts, which is assumed constant. . More pre- 
cisely o^X(t) [N - X(t)] is the number of contacts per unit time 
between affected and unaffected individuals. 

This model is more realistic than that of Bailey because it is true 
that in the case of a disease epidemic, the infectiousness of the path- 
ogenic virus may increase or decrease with the time of the process and 
the time since it infects the body of the individual . 

The equations (1.1) and (1,2) are difficult to solve in the general 
case, therefore the investigators made some assumptions and solve them 
for the following cases: 

a. The probability p(t,x) is a function of t alone 

b. The probability p(t,r) equals e _ ^ c 

c. The probability p(t,x) is function of X. only 

d. The probability p(t) is constant in a finite interval 

e. There is a lag in the transmission of contagion. 

Landahl [3] treated the spread of some information as a flow of 
'particles' which execute random motions over a population of individuals 
and which may multiply or disappear. Equations are derived for the 
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density of these particles, and the number of individuals through 
which the 'particles' have passed, is calculated. The results are ap- 
plied to the following cases: 

a. Uniform spatial distribution with multiplication factor decrea- 
sing with time because of loss of interest. 

b. Multiplication factor is constant, but the rate of spread de- 
creases with multiple hearings. 

c. One dimensional region with a small starting region with or 
without an absorbing barrier. An absorbing barrier corresponds to the 
case where the individuals listen to the information but do not repeat 
it. 

d. Two-dimensional region with absorbing barrier. 

e. Continuous sources of information within a small region in one 
dimension. 

f. Uniform spatial distribution in which individuals do not 
respond to more than one hearing. 

Prihar [4] studied a mobile communication network where the trans- 
mitting and receiving stations are in motion. He used the classical 
Kinetic Theory of gases to compute the expected total number of vehi- 
cles contacted per unit time in two cases: 

a. The search for contact is continuous and the antenna is omni- 
directional . 

b. The search is intermittent and the antenna is directional. 

Finally, Mattei [5] assumed the probability of contact between any 

pair of stations is constant throughout the network. 

Prob (contact between any pair) = p = constant. 

Several formulas were derived in closed forms. 
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The mean number of contacted nodes at the 1st step is 



Z(l) = (n - l)p. 

The mean number of contacted nodes at the 2nd step is 
Z(2) = (n - l)q [1 - (1 - p 2 ) n " 2 ]. 

The marginal probability of the number of new nodes at each step is 



Prob[Z(k) = z(k)] 




z(l) = 0 z (2) = 0 



ylkl.2) 

z(Fa) = 0 



|(^(?j) [1 - q z tO)]z(l) [ q z(0)]y(l)| 

leg]) [i - q z(i, ] z(2) 




qz(k-l) j z (k) 




. k 

y(k) = n -)T z(i) q = 1-p 

i^0 



(1.3) 



where n = number of nodes in the network 
k = order of step 

C. THE PROPOSED MATHEMATICAL MODEL 

Although with the model of Mattei, one could derive many quantities 
of the network in closed and neat forms, the model is not realistic in 
real life as pointed out in the beginning of this paper. The probabil- 
ity of contact between any pair of stations is far from constant. It 
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depends on their distances to the center of attack, the size of the 
electronic equipments, the geographical positions. In the event of a 
war, the capital, the missile silos, the industrial centers are more 
vulnerable than a submarine under the sea or a communication satellite 
in orbit. 

In this paper, a model, in which the probability of contact between 
any pair of stations could be varied at will, could be assigned to each 
link, once all circumstances have been taken into account, is proposed 
and solved. The solution is obtained by simulation on the digital 
computer. The computer program allows the user to assign the probabil- 
ity of contact between all pairs of stations. 
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II. SIMULATION ON A DIGITAL COMPUTER 



The aim of this chapter is to obtain, by simulation on 
the digital computer using the Monte Carlo method, the de- 
sired quantities of a general random communication network, 
in which the link probabilities could be all distinct. 

A. THE MONTE CARLO METHOD 

The Monte Carlo method consists of solving various 
problems by means of the construction of some random process 
for each such problem, with the parameters of the process 
equal to the required quantities of the problem. These 
quantities are then determined by means of the observations 
of the random process and the computation of its statistical 
characteristics which are approximately equal to the re- 
quired parameters. 

For example, the required quantity x might be the math- 
ematical expectation E(f|) of a certain random . The Monte 
Carlo method for determining the approximate value of the 
quantity x consists of an N-fold sampling of the value of 
the variable Vj in a series of independent tests, (that is, 
<1M 2, ... ) and the computation of their mean value 

^ a Hi + *)*.*•• * 

IV 
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If the number of tests is large, one has 



) = x. 



It is appropriate to ask the question why one does not 
try to derive an analytical expression for the quantity x 
from the properties of the problem instead of designing 
a random process which imitates the reality. The answer 
is that many problems in the physical world are so complex 
that they do not have closed formula for the desired quan- 
tities. 

Suppose, for example, that the required quantity x 
is the definite integral of a function f(y) taken over 
the interval (a, b) , that is 




Fig. 2.1. Illustrating equation (2.2). 
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If the expression of f (y) is not simple (or even worse, 
if the function f(y) does not have a closed formula but 
is known only point by point) , x could not be evaluated 
analytically. The Monte Carlo method may be used as will 
be described next. 

Suppose we have a device called a random number gen- 
erator which generates two random variables ^ and V) 
uniformly distributed ever the intervals (a, b) and (o, c) , 
respectively. The Compound probability density function 
is given by 



The pair (^,ty) will be generated N times and each time, the 
condition 



will be tested. If condition (2.1) holds, the point 
(^ i , vj i) in Fig. 2.1 is inside or on the boundary of the 





otherwise 



f(?i ) .< 



( 2 . 1 ) 



computed area. Let N* be the number of times when this 



condition holds in N tests, we have 




f(y)dy = c(b-a) A*” ) 



N-> no N 



( 2 . 2 ) 
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In order to obtain good results, the number of tests 
must be large, say, in the range between 1,000 and 100,000. 
This is only possible with the aid of digital computers. 

The problem of a random communication network where 
the probability of survival of each link is assumed constant 
has been solved adequately by Mattei [5]. If one changes 
this assumption and lets the said probability vary to fit 
certain conditions in reality, the problem becomes untrac- 
table, then the only feasible solution is to use a 
numerical method such as the Monte Carlo technique. The 
required quantities are: the average number of contracted 

nodes, the average number of newly contacted nodes at each 
step, the terminal reliability, the probability distribution 
of new nodes at each step. The method for obtaining these 
quantities will be described in the following section. 

B. THE FLOWCHART OF THE SIMULATION PROCESS 

Suppose that one has a network consisted of n stations 
connected by a certain number of links, it is a determinis- 
tic network in the sense that, if a message is sent by a 
station to adjacent stations which relay it to the others, 
one knows exactly how many and what stations have received 
it. Suppose a disaster, which could be a nuclear attack, 
an earthquake, etc., happens to the network. Before any 
investigation of the damage could be made, one does not 
know which links have survived. If a message is sent, it 
is difficult to acertain how many and which stations will 
get the message. The best one could do is to assign to 
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each link a certain probability of survival. Let be 

the probability of communication between station i and sta- 
tion j, the probabilities between all pairs of stations are 



described 


by an 


n x n 


matrix as 




r 0 


P 11 


P 12 


P v 
In 




P 21 


0 


P 23 


• * * * P 2n 


1*0 

II 












/nl 


P n2 




.Pn, n-1 0^ 



where P is called the matrix of probabilities. 

Since one is not interested in the case where a 
station sends back to iteself a received 
message, all diagonal elements of JP are zero. From a set 
of n stations, one could form 2 n ^ n ~^. networks. Fig. 2.2 
shows two of them for n = 5. 




Fig. 2.2. Samples of random network 



© 
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If the number of stations increases, say n = 21 stations, 
the population of possible networks is 2^ 2 ^ X 2 ^ 2 . 7 x 

1 y ft 

10 A , which is large. Suppose that one investigates the 

propagation of message via the originator which is the 
station 1. Starting with this station, one generates a 
random variable r uniformly distributed between 0 and 1, 
and check the condition 




step 0 step 1 step 2 step 3 

z (0) =1 z (1) =2 z (2) =4 z (3) =2 z(4)=0 



Fig. 2.3. A tree of propagation 



r 




(2.3) 



If this condition holds, there is a contact between station 
1 and station 2 and one increases the random variable z(l), 
which represents the number of newly contacted nodes at 
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step 1, by one unit. One generates again a random number r, 
tests condition (2.3) for P and increases z(l) by one 
unit if condition (2.3) holds, and so on until one arrives 
at the last station, which is station 21. One repeats the 
same procedure with each station of z(l) relative to the 
remaining network, the number of stations contacted at this 
second step is designated z(2). One repeats the same pro- 
cedure with each station of z(2). The procedure stops at 
the step when one could not contact any more stations. 

As an example, consider the network in Fig 2.3, stations 
5 and 17 have been contacted at step 1 so that z(l) = 2. 

The stations which have received the message are discarded 
from the network, these are stations 1, 5, 17. Next, to 
each station contacted at step 1, the same procedure made 
for station 1, is repeated. In the example, stations 2, 6, 

8 and station 14 have been contacted respectively by station 
5 and station 17 at step 2, z(2) = 4. Stations 4 and 13 
have been contacted by station 6 at step 3, z(3) = 2. The 
propagation of the message is over at step 4 because stations 
4 and 13 cannot carry the message any further. 

The tree in Fig. 2.3 summarizes the process in various 
steps. At any time the computer needs to keep the infor- 
mation of only one tree, which saves a lot of memory 
storage. To start the next sampling, one just has to erase 
all information of the first sampling run except, of course, 
the results z(l), z(2), .... The general procedure for 

the second sampling is the same as for the first. 
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After K samplings, the average number of new nodes at 
each step is given by 



z Ci) 

Z(2) 



Zj (i) + z 2 ( i ) + z 3 (i) + . . . + z k c i ) 

K 

Zj (2) + Z 2 C 2 ) + Z 3 (2) + . . . + Z k (2) 

K 



Z (KMAX) + Z (KMAX) + . . . + Z, (KMAX) 

Z(KMAX) = — = jr-= 

where KMAX = the highest order of steps in the tree. 

The average number of contacted nodes is 

X = 1(1) + 1(2) + . . . + Z(KMAX). 

Z 1 (l) + Z 1 (2) + . . . + (KMAX) 

+ Z 2 (l) + Z 2 (2) + . . . + Z 2 (KMAX) 



+ Z k ( 1 ) + Z k (2) + . . . + Z R (KMAX) 

- 

X(l) + X(2) + . . . + Y ( KMAX) 
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where X(i) is the average number of contacted nodes up to 
step i. 

The Kth terminal reliability is commonly defined as 



r(K) 



_ Z(K) 
n-1 



and the weak connectivity is defined as 



r 



X(KMAX) 

n 



Definition . In this paper the link density of the network 
is defined as the ratio of the number of links between sta- 
tions to the number of emitting stations. 

During the process of simulation, the probability dis- 
tribution of new nodes at each step is also calculated and 
displayed in an n x n matrix whose rows represent the num- 
ber of steps and columns represent the number of new nodes 
as follows: 



# of contacted stations 
0 1 2 3 ... n 

0 

1 

2 

. order of step 




Fig. 2.4. Probability distribution of the number of con- 
tacted stations at each step. 
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A flowchart of the simulation is drawn in detail in 



Fig. 2. 
program 

IRMAIN 



NEW 



ICOPY 



ZVEC 

XVEC 

GVEC 



The list of the variables used in the computer 
is described as follows. 



vector whose components 
If the ith component is 
has yet been contacted, 
ith has been contacted. 

vector whose components 
the ith component is 1, 
has just been contacted 
step . 



have the value 1 or 0. 
1, it means station i 
0 means the station 

r 

have values 1 or 0. If 
it means the station i 
during the preceeding 



vector with 1 or 0 as components, if the ith 
component is 1, it means the station i has been 
contacted during the current step. Before to 
go to the next step, one has to make; NEW = 
ICOPY. 



vector whose ith component represents the num- 
ber of newly contacted stations at step (i - 1) . 

vector whose ith component represents the aver- 
age number of contacted nodes at step (i - 1). 

vector whose ith component represents the 
(i - l)th terminal reliability. 
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RELIAB = 


vector whose ith component represents the 
(i - 1) terminal reliability. 

vector whose ith component represents the prob- 
ability distribution of new nodes. 


RNGE 


vector whose components determine the x, y 
ranges of the plot in the subroutine UTPLOT to 
plot the average number of new nodes. 


PRDIST = 


square matrix whose (i, j) component represents 
the probability that the number of contacted 
nodes is j at step (i - 1) . 


PROB 


square matrix whose (i, j) component represents 
the probability of communication from station 
i to station j . 


RLIAB 


square matrix whose ith row represents the ter- 
minal reliability for the ith case (different 
probabilities, different standard deviations, 
etc . . . ) . 


XXVEC 

\ 


rectangular matrix whose ith row represents the 
average message propagation for the ith case 
(different probabilities, different standard 
deviations, etc...). 
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KSTEP 

KHECK 

KHECK1 

SUMND 

ZNODES 

KMAX 

KKMAX 

CONECT 

NDEMIT 

LABEL 

MC 

IX 



integer which represents the order of the step. 

integer whose value is 1 or 0 used to initialize 

ICOPY to (0) or to compute the number of emitting 
stations . 

number of contacted nodes up to step K. 

number of contacted nodes during step K. 

highest order of step attained during the pro- 
cess of simulation. 

KSTEP + 2 used in the computation of PRDIST 
minor variable. 

represents the weak connectivity of the network. 

the number of emitting stations during one sim- 
ulation . 

used in the subroutine DRAW. 

integer used in the subroutine DRAW. 

integer used in the subroutine OVFLOW and RANDOM. 
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A 


= the average number of contacted stations at 

the 1st step in the model of Mattei or defined 
as n/A, standard deviation in the general model. 


S 


= the standard deviation of destruction in the 

general model, defined as: S = n/A. 


X and Z 


= Normalized variables, defined as X = I/S and 

Z = J/S. 


K 


= the number of tests in the process of simulation. 



P and PX = outputs of the subroutine NDTR which are areas 
under the normal curve. 

The following subroutines of the NPGS computer center library 
are used in the program. 



UTPLOT 


= to plot the average number of new stations curve. 


DRAW 


= to plot the probability distribution of con- 

nected stations, the terminal reliability and 
the average number of contacted stations. 


OVFLOW , 


RANDOM = used to get a random variable (0,1). 
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NDTR = to obtain the area under the normal curve as 
probability of communication between two sta- 
tions. 

To investigate the general model where the destruction 
is assumed normally distributed around the center of attack, 
block(2) has to be replaced by block(l). 
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SIMULATION OF A RANDOM COMMUNICATION NETWORK 




IRMAIN = (11 


. . . 1) 


ZVEC = (10 . . . 


0) 


NEW = (00 . . 


. 0) 


XVEC = (10 . . . 


0) 


ICOPY = (00 . 


. . 0) 


GVEC = (0 . . . 


0) 


KSTEP = Oj 


KHECK = 


0, KHECK1 ' = 0 




SUMND = 1.0, 


ZNODES 


= 0.0, KMAX = 


0 


KKMAX = 0, 


CONECT = 


0.0, DNSITY = 


0. 


NDEMIT = 0, 


MC = 1, 


IX = 102530469 



A=S=Z=X=P=PX= 0.0 

RNGE = (20.0, 0.0, 21.0, 0.0), K = 38,000 



I CALL OVFLOW 




Fig. 2.5. The Flowchart of the simulation. 
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PLOT: TERMINAL 


r 






RELIABILITY AND 




L ^ 3 


j 


1 

nDflDATATT J 







© 

i, OO 



L 1 



L c-L t- 1 



-<r 



PLOT PROBABILITY 
DISTRIBUTION 
CURVE 



\ 



BLOCK (I) 



BLOCK (2) 



A <- A f 1.0 PROB (I,J) = 








S <- 21/A PROB (J,I) 








X <.- I/S PROB (I,J) = 0 






ArA+ 1.0 


CALL NDTR (I = J) 






PROB <- ( A ) 
N-r 


Z <.- J/S 






PRDIST <•- (0) 


CALL NDTR 






PROB (I,J) <- (PX - 0.5) (P- 0.5) 


. _ 




; 






» 

‘ - - - - - > - - - 




PROB 





}co 



(D 



M <•— 1 


m<;k 




M C M+ 1 





PRDIST (1,2) = M 
IRMAIN (1) = 0 
NEW (1) = 1 
ICOPY (1) = 1 






© 



© 



© 



Fig. 2.5. continued. 
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DNSITY ^-0.0 
CONECT 0 . 0 
ZVEC (0) 
RELIAB <- (0) 





Fig. 2.5. continued. 
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V 



V 





Fig. 2.5. continued. 
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I' 



• KMAX > 


KSTEP ) 


vF 


KMAX 


KSTEP 



I 



ZVEC (KSTEP + 1) ZVEC (KSTEP + 1) + ZNODES 
PRDIST (KSTEP + 1, ZNODES + 1.0)*- 
PRDIST (KSTEP + 1, ZNODES + 1.0) + 1.0 
ZNODES d.- 0.0 



* 




L 

PRDIST(I6, 1) <- 
PRDIST (16,1) + 1.0 



Fig. 2.5. continued 




r 



NDEMIT = 0 



J 



DNSITY 

CSUMND - 1.0) 
NDEMIT 

NDEMIT 0 



= DNSITY 



KSTEP 0 
KHECK 0 
SUMND 1 . 0 
IRMAIN (1) 
NEW (0) 
ICOPY (0) 



Fig. 2.5. continued. 
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NEW r- I COPY 
KSTEP KSTEP - 1 







| 




/ 

f 

L 


KMAX 


<. KSTEP 


" 1 l ~\ 
> 






i T 






KMAX 


KSTEP 





KHECK 4 - 0 

ZVEC (KSTEP + 1) = ZVEC (KSTEP + 1 ) + ZNODES 
XVEC (KSTEP + 1) = XVEC (KSTEP + 1) + SUMND 
ZNODES 0 

PRDIST (KSTEP + 1, ZNODES + 1 ) PRDIST(KSTEP+1 , ZNODES + 1) 

+ 1.0 



Fig. 2.5. continued 
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C. COMPUTER PROGRAM 



The computer program is written to analyze a network of 21 sta- 
tions, however this number could be increased at will. First one has 
initialize the matrix of probabilities of communication between sta- 
tions to desired values. The probability distribution function could 
be discrete or continuous. 

For example, in the case where the probability of communication 
between any pair of stations is constant, it suffices to write 

Table 2.1. Fortran statements for the model where p is constant 



DO 16 I = 1, N 

DO 15 J = 1, N 

IF (I. EQ. J) GO TO 14 
PROB (I,J) = A/(N - 1) 
GO TO 15 

14 PROB (I,J) = 0.0 

15 CONTINUE 

16 CONTINUE 



If one wishes to analyze a network where the destruction is normally 
distributed around the center of a nuclear attack, it suffices to 
replace table 2.1 by the following table 
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Table 2.2. Fortran statement for a normal model 



S = 21/A 

DO 16 I = 1,N 

X = I/S 

CALL NDTR (X, PX, D) 

DO 15 J = I, N 

IF (I. EQ. J) GO TO 14 

Z = J/S 

CALL NDTR (Z, P, D) 

PROB (I,J) = (PX - 0.5) . (P - 0.5) 
PROB (J,I) = PROB (I,J) 

14 PROB (I,J) = 0.0 

15 CONTINUE 

16 CONTINUE 



The computer program could be used to analyze networks of any size 
N by redimensioning the variable names which are vector arrays, matrix 
arrays to the size of the network and enter the data card according 
to the changes. The program prints out the link density and the weak 
connectivity of the network, as well as the curves of the average 
number of contacted stations, the probability distribution function of 
the number of contacted stations at each step and the terminal relia- 
bility of the network. 
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SUMND= SUMND+ L 
ICCPV ( J ) = L 



ZNC'DES=ZN0DES+1 
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OTTING OF THE PROBABILITY DISTRIBUTION OF THE NUMBER OF NEW 
DES AT EACH STEP 
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III. COMPARATIVE ANALYSIS 



For the purpose of comparison, the computer program has been used 
to investigate a random communication network having 21 stations 
whose probability of contact is assumed constant. Three cases have 
been considered where the probability is equal respectively to 0.05, 
0.10, and 0.15. The sampling size is calculated to have an error of 
0.01 and the results are compared with the theoretical values. 

A. THE CALCUATION OF THE SAMPLING SIZE 

Suppose that one station tries to send the message to all other 
stations. The expected number of contacted stations is given by 




(3.1) 



The variance is given by 



<T = (n - l)pq 



where 



n = the size of the network 



p = the probability of contact 



q = 1 - p. 



Equation (3,1) becomes 




in the case where n = 21 and p = 0.05. 
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Let Z(l) designate the random number of contacted stations at the 
first step in one sample. Since the population is infinite, the mean 
and the standard deviation of the sampling distribution of means are 
given by 



hur z o) = h 



and 



°"~Z (1) 



a~ 




where N = the sampling size. 

For large values of N (N 30) the sampling distribution of means 
is approximately a normal distribution with mean and standard de- 
viation • The accuracy of the approximation improves as N 

gets larger (Central limit theorem) . 

Suppose that one would like the value of Z(l) differ from the 
population mean by 0.01 with a confidence level of 95.45%, one has 
to satisfy the inequality 





1 Z(1 )• -p- | 


2 °z(i) 4 °- 01 


(3.2) 


or 


2 \f(n-l)pq 
\[W 


^ 10 -2 




or 


N ^ 4.10 4 


. (n-l)pq = 4.10 4 . 


(21-1) (0.05) (0.95) 


then 


N y 3.8 x 


10 4 





The program has been run through 3.8 x 10 4 iterations. The results 
of the simulation are compared with the theoretical values in the fol- 
lowing section. 
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B. COMPARISON BETWEEN SAMPLING AND THEORETICAL VALUES 



a. Average Numbers of Contacted Stations at the First and Second 
Steps 

The average numbers of contacted stations are given by [5] 
Z(l) = (n-l)p 

Z(2) = (n l)q|l - (l-p 2 ) n " 2 }. 



They are calculated for three values of the probability p and compared 
with the results of the sampling in tables 3.1 and 3.2. The difference 
between sampling and theoretical results is small. 



Table 3.1. Values of Z(l). 



Probability 


Theoretical 


Sampling 


Difference 


0.05 


1 


0.9978 


0.0022 


0.10 


2 


2.0007 


-0.0007 


0.15 


3 


2.9870 


0.0130 

1 



Table 3.2. Values of Z(2). 



Probability 


Theoretical 


Sampling 


Difference 


0.05 


0.8825 


0.8708 


0.0117 


0.10 


\ 

3.1289 


3.1334 


0.0045 


0.15 


5.9676 


5.9370 


0.0306 
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The sampling size has been calculated for the value of probability 
equal to 0.05. Table 3.1 shows that the difference between the sampling 
and theoretical values Z(l) - 1 is 0.0022, less than 0.01 as it is 
expected. 

b. Average Total Numbers of Contacted Stations 

Fig. 3.1 gives the numbers of stations which will be eventually 
contacted for different values of probability. They are tabulated in 
Table 3.3 along with the corresponding theoretical values. 



Table 3.3. Average total number of contacted stations 



Probability 


Theoretical 


1 

Sampling 


Difference 


0.05 


4.50 


4.40 


0.10 

f 


0.10 


14.00 


14.07 


-0.07 


0.15 


19.00 


19.10 


0.10 



The theoretical [5] and sampling curves giving the average numbers of 
contacted stations up to step K are similar. Both show that the 
propagation of message is over between steps 4 and 6. 
c. Terminal Reliability 

The terminal reliability curves in [5] and in Fig. 3.2 are 
very similar. They show that the highest reliability occurs between 
steps 2 and 3 and is negligible from step 6. The terminal reliabil- 
ity is the measure of the speed of propagation of the message. 
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Table 3.4. Values of highest terminal reliability 



Probability 


Theoretical 


T 

Sampling 


1 

1 

Difference 

i — - 


0.10 


0.175 


0.175 


0.0 


0.15 


0.34 


0.32 

i 


0.02 



d. Probability Distribution of the Number of New Stations 
Contacted at Each Step 

Fig. 3.3, 3.4, and 3.5 give the probability distribution of 
new stations contacted at each step, they correspond to the formula 
(1.3), they are very similar to the theoretical curves in [5], 

In summary, the computer program has yielded results almost 
identical to the theoretical values if the sampling size is large, 
between 4 x 10^ and 10^. Since the program was written without any 
constraint, any assumption on the probability of communication between 
pairs of stations, it could be used with confidence to study any ran- 
dom network on the condition that the matrix of probabilities is 
initialized to values which reflect each particular circumstance. As 
an example, in the next chapter, the program is used to investigate 
a communication network whose destruction is supposed normally dis- 
tributed around a nuclear explosion or an earthquake. 



45 



a 
. * 
CJ 




Fig. 3.1. Average message propagation 

x scale = order of step (2 units per inch) 

y scale = number of contacted stations (3 units per inch) 
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Fig. 3.2. Terminal Reliability 
x scale = order of step (2 units per inch) 
y scale = terminal reliability (0.1 unit per inch) 
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Fig. 3.3. Probability distribution of the number of contacted 
stations at each step. 

Probability = 0.05 

x scale = number of contacted stations (3 units per inch) 
y scale = probability distribution (0.1 unit per inch) 
k = order of step 
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Fig. 3.4. Probability distribution of the number of contacted 
stations at each step. 

Probability =0.1 

x scale = number of contacted stations (3 units per inch) 
y scale = probability distribution (0.1 unit per inch) 
k = order of step 
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Fig. 3.5. Probability distribution of the numbers of contacted 
stations at each step. 

Probability = 0.15 

x scale = number of contacted stations (3 units per inch) 
y scale = probability distribution (0.1 unit per inch) 
k = order of step 
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IV. INVESTIGATION OF A GENERAL MODEL 



A. REPRESENTATION OF THE MODEL 



Consider a communication network of n stations which may be rep- 
resented by an n x n matrix whose elements are 1 or 0 as 



P, 



12 3 



0 10 0 1 

0 0 10 0 

1 0 0 0 . 
0 . . . . 

1 . . . . 



n 



n 

1 

0 

1 



(4.1) 



Suppose there is a perturbation near station 1, whose effect is to 
destroy the links of communication according to a normal distribution 
N(0,<r) around a center. This means that the destruction is maximum 
at the point of impact and decreases radially according to the 
equation 



P(x) 




_J 

cr yz-rr 




dt . 
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/ 




Fig. 4.1. Illustrating the equation (4.2) 



where p(x) = probability of survival of a link of a station 



x = distance of a station to the impact 



C~ = standard deviation of the destruction 



The probability of communication between stations i and j is 



Pij = P(Xi) p(Xj) 
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H 









1 



(4.2) 




<rY2TT 



t2 



*2 cr 4 




_t£ 

_J_’ e 20 2 

c-Orr 



dt. 



Suppose that the stations are numbered according to increasing dis- 
tance and the network is complete, then all elements of the matrix^ 
have to be replaced by the values of equation (4.2). The value of 
the standard deviation measures the extent of the destruction. The 
simulation was made for three cases where station 21 is one, two and 
three standard deviations away from the impact. 



first case 


cr = 21 


extensive destruction 


second case 


O- =fi= 10.5 


medium destruction 


third case 




minor destruction 



The subscripts of the stations which represent their distance to the 
impact, have been expressed in standard units and the subroutine NDTR 
has been used to initialize the matrix^. The sampling size has been 
taken equal to 3.8 x 10^. The results are summarized in the following 
section. 

B. RESULTS OF THE SIMULATION 

When the farthest station is one standard deviation from the im- 
pact, i.e, the case when the destruction is extensive, the message is 
not carried to many stations. On the average, 1.26 stations receive 
the message, the link probabilities being very low. In the case of 
medium destruction, the farthest station two standard deviations from 
the impact, 4.17 stations receive the message. The propagation 
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ceases at step 5. 7.45 stations receive the message when station 21 

is three standard deviations away from the impact. The above results 
are obtained with the assumption the communication network forms a 
complete graph, i.e., there is a link between any pair of stations. 

It is usually not the case in real life, therefore the above results 
are only the upper bound of number of stations which will receive the 
message in case of a war. 

If one compares with the model where the link probability is con- 
stant, one finds the latter model always gives much higher results 
than the current model which is nearer to reality. Table 4.1 summar- 
izes the results of two models. 



Table 4.1. Number of contacted stations in two models 
compared . 



Probability 


Mattei 


Standard Deviation 


Normal Model 


0.05 

1 


4.40 


21 ; 


1.25 


0.10 


14.07 


i 

10.5 1 

i 


4.17 


0.15 


19.10 


t 

7 

i 


7.45 



Fig.' 4. 2 shows the average message propagation in the three cases. 

The weak connectivity of the normal model is also smaller than the 
Mattei model, it is equal to 0.06, 0.20 and 0.35 when the farthest 
station is away from the impact one, two and three standard deviations 
respectively. The link densities are respectively 0.08, 0.36 and 
0.67. All results point out that in the normal model, the propagation 
stops rapidly. The terminal reliability in Fig. 4.3 is a monotoni- 
cally decreasing function in the first case, attains the maximum 
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value at step 3 for the second and third cases. It means that the 
propagation increases, attains the culmination then decreases and 
stops in average at step 5. The curves of probability distribution 
function in Figs. 4.4, 4.5 and 4.6 of the number of contacted stations 
at each step show also small values compared with the Mattei model. 
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Fig. 4.2. Average message propagation 

x scale = order of step (2 units per inch) 

y scale = number of contacted stations (3 units per inch) 
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Fig. 4.3. Terminal reliability 

x scale = order of step (1 unit per inch) 

y scale = terminal reliability (0.1 unit per inch) 
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Fig. 4.4. Probability distribution of the number of contacted 
stations at each step 

Standard deviation =21 

x scale = number of contacted stations (3 units per inch) 
y scale = probability distribution (0.2 unit per inch) 
k = order of step 
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Fig. 4.5. Probability distribution of the number of contacted 
stations at each step 

Standard deviation = 10.5 

x scale = number of contacted stations (3 units per inch) 
y scale = probability distribution (0.2 unit per inch) 
k = order of step 
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Fig. 4.6. Probability distribution of the number of contacted 
stations at each step 

Standard deviation = 7 

x scale = number of contacted stations (3 units per inch) 
y scale = probability distribution (0.2 unit per inch) 
k = order of step 
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V. CONCLUSION 



The analysis of the message propagation through a random communi- 
cation network has been studied by various investigators. To be able 
to solve the problem, each has adopted different assumptions on the 
probability of communication between stations. The main effort of 
this paper has been to solve the same problem without having to make 
any simplifying assumption on the link probability. The latter could 
vary to fit a particular circumstance. 

A flowchart has been drawn and a program for its implementation 
on the digital computer has been written. The program gives the de- 
sired quantities of a random communication network whose link 
probabilities could be varied to fit a particular situation by the 
Monte Carlo technique. The program has been tested with the model 
whose link probabilities were assumed constant and whose characteris- 
tics could be derived analytically. The results were found in 
agreement with the analytical results. The program has also been 
used on a model where the destruction was supposed to be normally 
distributed around the point of attack. The results also agree with 
the intuitive idea that, as the effectiveness of the destruction of 
the network increases, the propagation of message stops rapidly. 

In the computer program, the time taken by the message to travel 
between any pair of stations is assumed constant, which is not true 
in real life. Some stations take longer time to process and send the 
information than the others. It is proposed for future work that the 
time tjj , time taken by station i to process information and send it 
to station j, could be varied at will. With this second 
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generalization added, the simulation will give the average time of the 
propagation of the message and what routing procedure will give the 
fastest communication between various stations. 
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